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Abstract. Electromagnetically induced transparency (EIT) has been often 
proposed for generating nonlinear optical effects at the single photon level; in 
particular, as a means to effect a quantum non-demolition measurement of a 
single photon field. Previous treatments have usually considered homogeneously 
broadened samples, but realisations in any medium will have to contend with 
inhomogeneous broadening. Here we reappraise an earlier scheme [Munro et al. 
Phys. Rev. A 71, 033819 (2005)] with respect to inhomogeneities and show an 
alternative mode of operation that is preferred in an inhomogeneous environment. 
We further show the implications of these results on a potential implementation 
in diamond containing nitrogen-vacancy colour centres. Our modelling shows 
that single mode waveguide structures of length 200 jim in single-crystal diamond 
containing a dilute ensemble of NV — of only 200 centres are sufficient for quantum 
non-demolition measurements using EIT-based weak nonlinear interactions. 



PACS numbers: 42.50.Gy,42.50.-p, 42.50.Ex 

1. Introduction 

The importance of quantum mechanics to modern technology is indisputable. 
However, what remains to complete the 'quantum revolution' p] is the exploitation 
of coherent quantum mechanics in technological devices as well as the incoherent 
quantum mechanics responsible for, e.g., transistor electronics. Systems of strongly 
interacting photons and atoms have long been convenient systems for probing coherent 
quantum mechanics through the held of quantum optics. 

Of all the effects between coherently prepared atoms and light, electromagneti- 
cally induced transparency (EIT) is often promoted as an important building block 
for physics and device applications, because EIT allows the possibility of large optical 
nonlinearities accompanied by complete transparency [2]. EIT is a coherent quantum 
phenomenon whereby the absorptive and dispersive properties of a three (or more) 
state system can be tailored by using applied electromagnetic fields, and we discuss 
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its properties more fully below. First observed by Boiler, Imamoglu and Harris [3J, 
some of the proposed applications for EIT include magnetometry |1] , high-efficiency 
UV generation photonic switches [7], and optical [5] quantum gates, and light 

storage in first-in first out (FIFO) networks [HI OH]- Although the medium of choice 
is usually a vapour cell (e.g. Rb [11]), future technology may be more easily realised 
with solid-state media. EIT has also been studied in solids [T2] . fT3l [Ml \15\ [TB I \T7 \ [18] . 
magneto-optical traps [19] and Bose-Einstein Condensates |20j . 

Here, we concentrate on the possibility of using the lossless Kerr nonlinearity 
assocatiated with EIT for realising a quantum non-demolition (QND) measurement. 
QND via the cross-Kcrr effect between two distinct optical modes was originally 
proposed by Imoto, Haus and Yamamoto [21], and invoking the EIT induced Giant 
Kerr nonlinearity is a popular suggestion for realising a tt/2 phase shift for such a 
measurement. The idea that such a QND gate could be realised by weak nonlincarities, 
such as are routinely found in EIT systems, without a full n/2 phase shift induced on 
the detection beam, was introduced in Ref. [55] . In this system, the weak nonlinearity 
is effectively enhanced by the presence of a strong probe beam. However the earlier 
proposal did not consider all of the limitations of realistic systems, and in particular 
did not consider the effect of inhomogeneous broadening of the EIT medium: we do 
so here. Although we are concentrating on QND measurements, it has been shown 
that QND measurements effected by weak nonlinearities can act as a primitive for 
other quantum gates, [23] and this directly leads to the Qubus [24] [25] [26] [27] and 
related schemes for quantum repeaters [28] and cluster-state generation [29]. The 
results presented here should be equally applicable to these, and other EIT-related 
schemes. 

EIT is a well-known mechanism for generating optical nonlinearities without loss 
[2], for a recent review see Ref. [32 . The typical (and minimal) system in which to 
observe EIT is a three-state system in the A configuration with two driving fields, 
depicted in Fig. QJa) with states labelled |a), \b) and |c). EIT is a manifestation of 
quantum interference in an atomic system: Considering field 1 as a 'pump' and field 
2 as a 'probe', absorption of the probe is suppressed due to the coherence induced by 
the pump. The coherence gives rise to interference between the dressed states on the 
\a) — \b) transition, and hence a dip in the probe absorption. Because the absorption 
dip is due to a coherent two-photon resonance (the resonant two-photon transition 
from \a) to |c)) it can, in principle, be much less than the optical linewidth of the 
\b) — \c) transition, limited only by the decoherence on the \a) — |c) transition. 

Although transparency is the eponymous feature of EIT, of more practical interest 
is the steep and controllable dispersion curve that is associated with the transparency 
point. A typical example of the dispersive and absorptive spectra associated with 
EIT are shown in Fig. [2] At the transparency point (two-photon resonance between 
the ground states) there is a linear dispersion. This feature has led to dramatic 
demonstrations of ultra-slow group velocity light [5D] [M] [35] and is the basis for the 
Giant Kerr nonlinearity [33] in the N configuration of Fig. [T] (b). Considering the 
N system as a A system perturbed by an off-resonant transition, one can see that 
the detuned off-resonant third field induces a small light shift to |c). Although this 
is a small effect, because of the steepness of the dispersion, a large effect on the EIT 
resonance is observed. This property was proposed for achieving photonic blockade 
in cavity QED systems [35], and has undergone extensive theoretical investigation 
(e.g. Refs. [33 [SHIES SOI [41]) and recently observed [42] (although in the two-state, 
rather than four-state configuration). When field 3 is resonant, an absorptive, rather 
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Figure 1. (a) Schematic of the three-state A system under consideration. States 
| a) and |c) are ground (meta-stable) states, and have no direct spontaneous 
emission pathways. \b) is an excited state, which decays to \a) and |c) with equal 
probabilities. The \a) — \b) transition is driven by field 1, with detuning A and 
Rabi frequency Oi (pump). The \b) — \c) transition is driven by field 2 (probe) 
with detuning A + <5 (i.e. <5 is the shift from the mutual detuning A) and Rabi 
frequency O2 ■ (b) Schematic of four-state system in the N configuration, where 
the A system is modified by an interaction with field 3 driving (off-resonantly) 
the |c) — \d) transition with Rabi frequency C3 and detuning A3. To first order, 
we expect the |c) — \d) transition to perturb the EIT in the [a) — \b) — |c) system 
via the usual light shift, which plays a role equivalent to the two-photon detuning 
in the three-state scheme, (c) Proposed configurations for realising the QND 
measurement. Fields 1 and 3 enter a single-mode diamond waveguide constructed 
from a photonic bandgap material with tailored group velocity for field 3 and of 
order 100 jim long. Two illustrations of potential structures are shown, the lower 
follows Krauss 1301 , and the upper is via coupled cavity arrays in the style of Altug 
and Vuckovic 1311 Field 2 (not shown) illuminates from the side. The unknown 
field, field 3, propagates unperturbed, but field 1 is phase shifted by the number of 
photons in field 3, and this phase shift can be read out using heterodyne methods 
(not shown). 
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Figure 2. Typical dispersion (a) and absorption (b) curves associated with EIT 
in a A system (solid lines) compared to those of a two-state system with equal 
probe couplings (dashed lines). The narrow EIT resonance is associated with 
linear dispersion at 8 = 0. The steep slope makes it easy to perturb the EIT 
resonance and gives rise to the Giant Kerr effect. 



than dispersive, nonlinearity is observed, which has been studied theoretically [7] and 
experimentally [331 S3] , and will not be treated here. Neither do we consider the many 
related atomic configurations, e.g. the Tripod [45], extended N [46], Chain-A [47], and 
M-scheme [3H], although they all offer potential improvements over the conventional 
./V scheme. 

One important detail for realising nonlinear interactions in the N system is that 
of group velocity matching. In a travelling wave geometry to realise the N system 
of Fig. Q] (b), we require the pulses that describe fields 1 and 3 to be temporally 
coincident for the maximum cross-Kerr interaction (the classical pump field 2 can 
be assumed to be derived from a large uniform field, and so is exempt from this 
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criterion). Because field 1 is travelling under EIT conditions, it will be propagating 
with extremely slow group velocity: field 3 is not. We see that the amount of mutual 
interaction would therefore be expected to be limited by the temporal walk-off of 
the two pulses. There have been many suggestions in the literature to counter this 
effect (e.g. Refs. [49l [50]) which invoke varying levels of complexity of the interaction 
medium. It is also possible to control group velocity by modification of the medium, 
e.g. by using tailored photonic bandgap structures [5TJ[52]. Because we are mainly 
interested in solid-state implementations, we will assume that the system is embedded 
in a photonic-crystal structure where the group velocity as seen by field 3 is tuned 
by the structure to balance the EIT induced group velocity seen by field 1, and we 
therefore will not treat this detail further. 

Although much of our treatment in this paper will be system independent, it 
is important to note that a major motivation for performing this reappraisal of weak 
nonlinear gates is the availability of a new material for observing optical EIT: diamond 
containing the negatively charged Nitrogen- Vacancy colour centre (NV). This material 
has shown quite remarkable results, including single photon generation (e.g. [53 ), 
room-temperature Rabi oscillations [54] , and spin-spin coupling [55] [56] [57] [58] . EIT 
has been demonstrated in NV diamond in the rf [13] and optical regimes [16[ [17] . 
The intrinsic properties of NV centres can also be modified by fabricating optical 
structures directly in ultra- nanocrystalline [33] or single crystal diamond [5D] [ST] [55] , 
by growth on preexisting optical structures [63 , and also by Stark shifting which 
has been demonstrated on bulk samples [64] and spectrally resolved centres [65] [66] . 
The time is therefore ripe to examine NV diamond for the goal of optical quantum 
information processing. 

In the next section, we discuss EIT and the properties of coherently driven A 
systems. By including inhomogeneous broadening and expanding about the EIT point 
we are able to derive analytical results for the absorption and dispersion without 
the usual assumption of weak probe and strong pump fields. We are also able to 
determine optimal ratios for pump and probe Rabi frequencies. In Section [3] we show 
the results for cross-Kerr nonlinearities in the four-state N system, and in Section 2] 
we present our main results, which use the results from the preceeding sections to 
design structures based on diamond containing the Nitrogen- Vacancy colour centre 
that should be sufficient to realise a number-discriminating Quantum Non-Demolition 
operation. 

2. Three-state A system 

A schematic of our model three-state system is shown in Fig. [T] Following Shore [6TJ , 
we write down the Hamiltonian under the rotating wave approximation in matrix form 
with state ordering \a), \b), |c) 

H = h (Aa bb + 5a cc + Qia ba + £l 2 o- c b + h.c.) , (1) 

where field 1 (field 2) drives the \a) — \h) (\b) — |c)) transition with Rabi frequency f2i 
(f2 2 ) and detuning Ai (A 2 ), and er^ = K)0'|- As we are operating near the two-photon 
resonance, we set Ai = A and A 2 = A + S where A is the mutual detuning, and <5 
is the detuning from two-photon resonance. Spontaneous emission at rate T is from 
1 6) to | a) and |c) with equal probabilities. We treat inhomogeneous broadening by 
considering a distribution of A, i.e. the mutual detunings. As we are only considering 
inhomogeneity on the excited state distribution, the two-photon detuning, 5, does not 
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vary. Note that although we will use terminology such as 'pump' and 'probe', we will 
usually make no assumption about the relative strength of these fields. In this way 
our analysis is analogous to the cases treated by Wielandy and Gaeta on EIT in the 
strong pump regime [68] or the parametric EIT regime [69j . 

One way to proceed in gaining insight on the EIT problem set out in Eq. [1] is to 
construct the master equation and determine the steady state solution. This can be 
expressed as 

3 

where we have introduced the usual density matrix, p and the Liouvillian super- 
operator, C[B,p] which describes the effect of the generalised decoherence channel B 
with rate T > on the density matrix, and is summed over all decoherence channels. 
The Liouvillian operators are defined 

C[B. p] = BpB^ — - (B^Bp + pB^B) . (3) 

In our case, we restrict ourselves to the case that the system is limited by spontaneous 
emission, and that the decoherence between the ground states can be neglected. This 
approximation is warranted because of the very long ground state decoherence rates 
in important systems of interest (e.g. diamond containing the negatively charged 
Nitrogen- Vacancy colour centre |56J, or rubidium vapour cells [34]) but does limit the 
minimum Rabi frequencies that can be applied to be larger than this decoherence 
rate. So here, the B will be the one-way (spontaneous emission) transitions from \b) 
to either \a) (cr a b) or |c) (a c b) with rate, T/2. 

To analyse Eq. [5] in the steady state, we convert the master equation into 
superoperator form, and write 

'p={-iV + C)p, (4) 

where pis the vector obtained by writing out the density matrix elements, and V and C 
are the superoperators describing Hamiltonian and decoherence processes respectively. 

In addition to the master equations, we also include the effect of inhomogeneous 
broadening. Inhomogeneous broadening has been treated previously in the context of 
Doppler-broadened EIT in vapour cells (e.g. [7(11 EI]), an d is treated as a Gaussian 
distribution of the absolute energy of \b), which in turn is manifested as a variation 
in A across the sample, i.e. there is a probability distribution of detunings 

(A- A„) 2 " 



2 7 2 



(•5) 



where A = J™ AP(A)dA is the mean mutual detuning with respect to the 
inhomogeneous linewidth, which has standard deviation 7. Fig. [3] shows spectra with 
the real and imaginary parts of the \b) — |c) coherence, pb c (proportional to the probe 
dispersion and absorption respectively) with £7i = T/2 and 0,2 = T/10 to illustrate the 
effect of inhomogeneous broadening on the resonant EIT profile. Figs. [3] (a) and (b) 
show the absorption and dispersion curves for A from to positive values showing the 
effect of increasing A. Note that spectra of the same colour go together. Although the 
overall spectra are quite different, in the vicinity 5 = they are all locally similar. This 
is clearer in Figs. [3] (c) and (d) which shows spectra from large negative to positive A 
superimposed, showing the local similarity strongly. Figs. [3] (e) and (f) show the effect 
over the whole inhomogeneous linewidth with 7 = 10r (solid blue line) compared with 
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Figure 3. EIT in for homogeneously broadened systems. (a) Re(p cb ) as 
function of <5 for the mutual detunings of A = 0, 2r, 4r (blue, black and red 
lines respectively) to illustrate the effect of moving away from resonance, (b) 
Absorption [Im(p c {,] under the same conditions. Again note that all of the EIT 
transparency windows overlap, (c) Closeup in the vicinity of 8 = 0, showing 
Re(p c ;,) this time for A = — 2r, — T, 0, T, 2r (blue, green, red, cyan and 
magenta lines respectively) (d) Analogous spectra for Im(p c (,). 



the a sample with equal total population but only homogeneously broadened sample 
(black dashed line). Note that these plots may equally well be interpreted as classic 
hole-burning spectra [72] • The self-similarity of the EIT traces is not dependent upon 
the mutual detuning, as is illustrated in Fig. [4] which compares Ao = with Ao = — T. 

The self-similarity of the EIT profiles about 6 = may be understood by 
considering pi, c in the steady state. Setting p — 0, we determine the null-space which 
gives the non-trivial steady state solution for /^ ss ' . In general the analytical results 
are somewhat complicated. By confining our interest to the region around 6 = 0, i.e. 
in the vicinity of the two-photon resonance, we may perform a series solution for p( ss ) 
in powers of 6, which yields to first order 



+ «i ~ (ft? + ft|) 



Po.a = o2 | o2 - — -Z^S, (6) 



^fb [a (n; - nl) - it (n; + nf)] * fa . 
p ---WTv 2 + ^T^f 5 - p - (8) 



Pbb = 0, (9) 



(ii) 
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Figure 4. Comparison of EIT traces as a function of the detuning from the 
two-photon resonance condition, S, for one-photon resonance, (a) real part, (b) 
imaginary part, and finite mutual detuning, A = 2Q, (c) real part and (d) 
imaginary part. In each case the dashed line is the single-atom result, whilst the 
solid lines are the average over the inhomogeneous line, with 7 = 10r. Although 
the overall lineshapes vary with A and inhomogeneity, the <5 = responses are 
equivalent to second order. 



Cursory inspection of these steady state results provides some very important 
properties of the EIT condition. Modelling of inhomogenous broadening was to be 
by varying A over the ensemble, however assuming field 1 is the pump and field 2 the 
probe, then the parameter of interest for EIT is p cb . From the above, we see the well- 
known linear dependance with detuning expected for an EIT resonance, which has no 
dependence on A, indicating that inhomogeneous broadening will affect neither the 
absorption nor dispersion seen by the probe field to first order in the probe detuning 
(but to all orders in the mutual detuning). This is an important result, highlighting 
that EIT is extremely robust to inhomogeneous broadening in the excited state. To 
see effects due to the inhomogenous broadening of the line, we will need to go to higher 
orders in 5. 

The total coherence is obtained by integrating over the inhomogeneous linewidth, 
so we have 



Qcb 



-f 

J — ( 



P(A)p cb olA, 



(12) 



where we have introduced g cb as the coherence integrated over the inhomogeneous 
line. As P(A)dA = 1, to first order in S, we have trivially that g cb = p cb . To see 
effects due to the inhomogeneous line, we must go to second order in 8, where we have 



o (2) - 

Pcb — 



(n 2 + n 2 ) 2 



-5 + 
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AS 2 



(n 2 + n 2 ) 2 
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(13) 
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Performing the integration over the inhomogeneous line yields 

/>OC 

/ P(A)pW dA , (14) 



r (n? - 30|) a j 2 r<5 2 

-o + - — 77^5 ^tT^ h 



(15) 



Regions with high dispersion are associated with regions of high nonlinearity, and 
hence will inform us in our search for optimal working points for our gates. The probe 
dispersion is defined 

R cb - (16) 
do 

and using the first order solution in eq. [TO] we obtain 

— (17) 

Note that if we set fii = \/3^2, then the second order correction to the refractive 
index in Eq. [13] is nulled, and we must therefore go to third order in S (or higher) 
to observe terms depending explicitly on the inhomogeneous broadening. This result 
also carries through to the dispersion, i.e. that the gradient is largest when Q± and 
f^2 are smallest and maximised when f2i = a/3^2- A graph showing log 10 (— R) as 
a function of Qi and f2 2 is shown in Fig. [5] along with a line showing the maximal 
dispersion at ill = \/3^2- These results may be useful to optimise slow and stopped 
light experiments, however for the purposes of QND measurements (as will be seen in 
Section [4| the susceptibility is simply maximised for minimum possible f2i. 

To explore the absorption we take the imaginary part of the solution for p c b, 
which is to third order in 5 

A = WTW " 2Ar (nt + n^r 6 + ° [S] - (18) 

Ignoring the third order correction, here we see the familiar quadratic dependance with 
respect to detuning, and as presaged above, there is no contribution to the absorption 
from A. We can infer the bandwidth of the EIT medium, by considering the effect 
of the EIT window on a pulse with finite bandwidth. Fig. [6]Ja) shows schematically a 
transform limited Gaussian pulse propagating through the quadratic EIT window. If 
we assume that the pulse is defined by some spectral width /, centered around the EIT 
window at S = 0, with the functional form exp[— <5 2 /(2/ 2 )]/v'27r/ 2 , then to deter mine 
the total (single atom) absorption we integrate the pulse over the EIT window, i.e. 

exp H 2 /(2.f 2 )] r; 2 f2 2 r 
(tti + ft 2 ) 3 
njn 2 rf 2 

~ (ft 2 + fi 2 )3' (20) 

and this quantity will prove essential in determining bandwidth requirements in the 
design of practical nonlinear gates. To explore this, we calculate in Fig. EJb) the 
bandwidth in units of the spontaneous emission, f /T as a function of Qi/T and f22/T 
that will achieve a A = 1%. Such analyses as these allow us to determine gate speeds 
for QND measurements and will be exploited in Section [3] 



-4= / 1 7=T J ,„/J,J > (19) 
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Figure 5. First order probe dispersion, plotted as log 10 (— R c b) as a function of 
Hi (pump) and Q2 (probe) in the vicinity of two-photon resonance in a three- 
state system. Note that the dispersion as seen by the probe will be largest when 
f!l = \/3f22, and goes as R c i, = — 3v3/(16fii) in this limit. This optimal solution 
for R is indicated by the solid line. 
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Figure 6. (a) Schematic showing the absorption seen by an optical pulse with 
finite bandwidth. To second order, the EIT window has absorption that is 
quadratic in 8, where a transform limited pulse is a Gaussian in 8, so the residual 
absorption per unit length seen by a pulse travelling through an EIT window will 
be found by taking the overlap integral of the pulse and absorption, (b) Log base 
10 of the maximum bandwidth in units of T that will attain A = 1 % as a function 
of f2i/r and f^/T. Note that although the bandwidths can be quite small, they 
are still of order the Rabi frequencies of the fields, which accords with our intuition 
about the width of the EIT window. The smallest allowable bandwidth coincides 
with the largest dispersions, again as one would expect. 

Finally, we comment on the group velocity seen under conditions of EIT. Group 
velocity reduction is one of the most dramatic consequences of EIT and is observable in 
the usual pump-probe arrangement (e.g. [201134] ). Perhaps surprisingly, large changes 
in the group velocity between fields (group velocity mismatch) can actually lead to 
a strong reduction in effective coupling. As mentioned previously, we propose group 
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velocity engineering as the solution to this mismatch, but it is essential to understand 
the group velocity under conditions of EIT in order to specify the propagation 
properties of the unknown field. 

The group velocity of a probe field can be determined using Eq. \T7\ so we have 
c 



9ab n -I- < >x "" 
'1^2 dA 



(21) 
(22) 



Xab — ^TTPab, (,2oj 



where rj is the bulk refractive index, and we have introduced the susceptibility 
2nAf nl 
eo£R.h £1 

d X ab _ 2llN \l\ 



and so 



OS £q£rTi 



'-Rab, (24) 



^^(nf + nl) 2 ' (25) 

where A/" is the number density of atoms. So in the limit that the group velocity 
reduction is large, we have 

ce Q e R h (n 2 + n 2 ) 2 



^ ■nujNn 2 ab n 2 

Ultimately we will be interested in the group velocity associated with quantised 
fields, rather than the semiclassical form above. Making the substitution 51 1 y / n7 = ^i 
where n\ is the number of photons in the mode, gives the group velocity seen by a 
mode with n\ photons of 

] 2 
L 2 , 



(26) 



ce e R H fo a m+ tf 

1TU)NlJL 2 ab Q 2 . 



(27) 

'2 

These results will be used in Section 01 especially with quantisation of the probe 
field. In general one would seek the largest possible v g that enables effective coupling, 
so as to minimise the reduction in the group velocity on the unknown signal which is 
not travelling under conditions of EIT reduced group velocity. Also the group velocity 
dispersion which will manifest with uncertainty in the number of photons in the probe 
field is also a potential source of error and should be minimsed. This implies a rule of 
thumb, that we should seek operation in the limit of large n\ to minimise the relative 
variation in v a . 



3. Four-state N system 

The four-state N scheme is one level structure that clearly shows a cross-Kerr effect, 
and is illustrated in Fig.[]Tb). There is much freedom to choose which fields correspond 
to pump, probe and driving, and all appear to have been treated in the literature in 
various places. For concreteness, we will treat the N system as a A perturbed by an 
off-resonant transition. Furthermore, we will be considering the effect of the |6) — \d) 
transition (hence field 3) on field 1, so in this section, the parameters of interest will 
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be p ab and g ab . For this case we set the operating point for the A system as <5 = 0. 
Under these conditions, the Hamiltonian is 

H = h (Aa bb + A 3 a dd + VL X (j ba + tt 2 v c b + ^o- dc + h.c.) , (28) 

where the (in general unknown) Rabi frequency of field 3 is £l 3 and it is detuned from 
the I b) — \d) transition by A 3 ^> fi 3 . 

In the case that the |c) — \d) transition is only homogenously broadened, we may 
treat the effect of field 3 on the \a) — \b) — |c) system quite simply. Our treatment here 
follows and extends Ref. [5]. The effect of a field 3 on the \c) — \d) transition can be 
seen as an off-resonant light-shift, which in turn perturbs the \a) — \b) — |c) A EIT. 
The strength of the light shift can be directly equated with the —5 (the shift will be 
in the opposite direction) from the previous analysis, i.e. 

A3 

in the limit that A3 3> f^. Recalling that the residual population in \d) is a source of 
error, we will assume that this population must be kept below some threshold, e, i.e. 

^=g) 2 <, (30) 
Using the previous result for the steady state of p ab , and substituting for <5, we have 

(Oj + SI2) A 3 

The Hamiltonian of Eq.[5Hlcan also be attacked using the superoperator approach. 
We assume that all decay from \d) is to |c) at rate T, and keeping the other terms the 
same from the A analysis, except A = 0. We can obtain the steady state response 
quite easily, but for clarity, we only report the coherence on the \a) — \b) transition, 
which is 

« = -(5A7W?W + 0[!y4 - <32) 

which is qualitatively very similar to the previous result, but with the inclusion of an 
extra absorption term. 

To move to the case of a finite inhomogeneous linewidth, we first recall that we 
may safely ignore linewidth to first order on the A system, however, we cannot do so 
on the \b) — \d) transition, as A3 appears directly in the coherence. Therefore, using 
the light-shifted treatment from Eq. [31] we write down the ensemble coherence as 



g ab = I P(A 3 )p ab dA 3 , (33) 
O x n|fil [°° exp[-(A 3 - A ) 2 /(27 2 )] 




dA 3 , (34) 



(35) 



where erfi(a;) = erf (ix)/i is the imaginary error function. Note that although we have 
treated the inhomogeneous linewidths for \d) as equivalent to that of \b), the analysis is 
practically unchanged if we took a more general case. It is instructive to examine the 
behaviour of just the exponential and imaginary error function terms, to determine the 



Single photon quantum non- demolition in the presence of inhomogeneous broadening. 12 




o.o — F — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — <- 

123456789 10 

d 



Figure 7. Graph showing optimal mean detuning from resonance in units of 
the inhomogeneous linewidth. Here we have taken our figure of merit to be 
J(d) = exp(— d 2 )erfi(d) which is the term in g a (, which depends explicitly on 
d = A/(2"/). Note that the results for d < 1 will not be relevant for our 
discussion because the condition f23 <C A3 may be unsatisfied, and there may 
also be significant resonant one-photon absorptions of field 3 



optimal ratio of mean detuning to inhomogeneous linewidth. Setting d = Aq/ \J2^ 2 
as the mean detuning in units of the linewidth, and J(d) = exp(— d 2 )erfi(c?) in Fig. [7] 
we plot J(d) vs d, as for constant 7, all the other terms are constant. The maximum 
of J can be seen as being just before d = 1, however, care must be taken in this limit, 
and there will still be appreciable absorption of field 3 in this limit (as there will be 
on-resonant atoms). However, detuning by 5 inhomogenous linewidths (i.e. d = 5) 
only has the effect of reducing the effective g a b by about one sixth compared to a 
homogeneously broadened sample at such a detuning: the larger effect is the penalty 
in having to go to such large detunings to avoid the line (i.e. for a homogenously 
broadened sample, one could work closer to resonance). 

An alternative transition to be considered for the readout is the |c) — \d) transition. 
The coherence associated with this can also be determined and is 



(2 a 3 + ir)(n? + of) 



Pcd = - 7ZX I ;r w ni I 7vF\ +0[n 3 ] J . (36) 



Integrating this over the inhomogeneity on the |c) — \d) transition does not yield 
analytic solutions, however by replacing the denominator by the large detuning 
approximation (i.e. the term (2A3 + iT)/2 is replaced by A3, we get 

Scd = - J^" 3 exp (-*L) erfi ( *L| . (37) 

For completeness, we also report the coherence on the |c) — \b) transition, which 

is 

2£l^£l2^3 

(2A 3 -»T)(n? + nif 

As the response here is clearly analogous to that of p a b we will not repeat any further 
discussion of using this transition for monitoring any probe state, except to note 
that changing between \a) — \b) and |c) — \b) may perhaps be useful for reasons of 
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experimental convenience in certain implementations, but otherwise would appear to 
hold no benefits. 

It is interesting to note that neither Eq.l3l?lnorl3l^exhibit self-Kerr effects, with the 
first self nonlinearities for p a \, appearing at fourth order, and for p cc i appearing at third 
order. The presence of Kerr terms would hamper state discrimination in nonlinear 
gates (see next section) [73] [74]. One should note that the canonical method for 
generating self-Kerr terms in EIT media is to allow a field to interact with more than 
one transition (e.g. [361 HS1 EZ] ) which effectively converts the cross-Kerr nonlinearity 
into a self-Kerr nonlinearity. The suppression of self-Kerr terms is desirable, as they 
give rise to pulse distortion which can limit the effectiveness of any gate based on 
nonlinear interactions. 

Finally, we comment further on the populations in the excited states. In Eq. [30] 
we presented a simple two state argument for the population in state \d), which was 
viewed as a potential source of error. We now examine the full solutions for p^ and 
Pdd, which turn out to be qualitatively similar to the analysis based on perturbing the 
EIT structure by the extra transition. 

Starting with p^d, we find that the expansion to third order in £1% yields 

te= O^A^OT + °^ (39) 

This result can be immediately interpreted as the usual, off-resonant population from 
Eq. [3D] (with spontaneous emission explicitly included) , scaled by a factor due to the 
diminished population in |c) because of the coherent population trapping in the A 
system. Similarly we can calculate the population in \b). Although the unperturbed 
EIT condition leads to no steady-state population in \b), the perturbed EIT will give 
rise to non-zero population. As above, this can be calculated and to third order in f^3 
we obtain 

In general, if we have SI2 = -\/zQ.i, then we will have 

z n 2 

Pbb = (z+l) 2 A| + (r/2)2' (41) 

1 n 2 

Pdd ~ z + l A| + (r/2) 2 - { ' 

Any population occupying the excited states will in general be available for 
decoherence and give rise to errors. By inspection of Eqs. [41] and [42] we observe 
qualitatively the same familiar criterion from Eq. 1301 (with minor corrections). 



4. Implications for the design of QND weak nonlinear detectors 

Our focus is on the construction of a device capable of achieving a QND measurement 
of the number of photons in a weak field. In this section we will combine the previous 
analyses above with realistic parameters that are achievable using a solid-state slow- 
light EIT waveguide, which contrasts with the more general discussion of EIT based 
nonlinear interactions. Our analysis focusses on QND measurement of a weak field 
with unknown photon number and we show discrimination between 0, 1 and 2 photons 
in the unknown field. In some of the parameter ranges discussed below, we also observe 
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distortions of the probe field. Whilst this is not a problem for QND discrimination, it 
may restrict the utility of our discriminator for use in quantum gates [25] . where the 
requirement for sequential use of the QND probe favours regimes where the Q-function 
of the probe is only rotated, and not also distorted by the nonlinear interaction. 

We will first describe the appropriate metrics for evaluating the performance 
of any such gate in a material independent fashion, and then conclude by presenting 
realistic operating conditions in several potential implementations as a guide for future 
demonstrations. For clarity, throughout this section we will restrict ourselves to the 
case that the \a) — \b) transition is probed by a weak coherent state, and the |c) — \d) 
transition has the unknown signal field. 

One of the most important parameters to determine the strength of the 
measurement signal is the effective Rabi frequency of the unknown pulse. Analysis 
of this will show that spatially confined structures (e.g. waveguides) will have 
a considerable advantage over free-space implementations. Following [75j we may 
express the single-photon Rabi frequency as 
Ax 



o 3 = -^vsrjjvi (43) 

47T 

where we have introduced A12, the resonant transition free-space wavelength, I the 
length of the medium (or with imperfect group velocity matching it will be the length 
of the effective interaction region) and recall that / is the bandwidth of the single- 
photon pulse. The Rabi frequency is then 

n 3 = fisV^, (44) 

and the other Rabi frequencies can be defined similarly. Note that there is usually 
a dependence on the beam waist in Eq. 1431 however this is compensated by the 
number of interacting four-state systems within the single photon spot size (so a 
larger spot interacts with more systems but with less strength). However, waveguide 
structures will still have significant advantages in minimising the beam cross-sectional 
area compared with free-space structures, with the effective medium length in free- 
space being ultimately limited by the Rayleigh range of the beam. 

We first need to connect the microscopic description presented above with the 
macroscopically observable quantities. In particular, when considering the four-state 
system with the probe field on the \i) — \j) transition, the phase shift seen by the 
probe [75] 

Kij = QiQij, (45) 
and the evolution of a state, \ai), impinging on the |i) — \ j) transition for time t is 
|c4) = exp (iKijt) \oi). (46) 

Before exploring numerical examples, it is instructive to study the nonlinear 
optical processes. In particular, it is natural to consider two transitions to be probed 
to effect quantum non-demolition measurements, i.e. we could probe the \a) — \b) 
transition, with the unknown field on \c} — \d) transition, or probe the |c) — \d) transition 
with unknown field on the \a) — \b) transition. Substituting for the density matrix 
elements (and assuming is a classical pump field) we have 

Kai ..^ masm^, j(^\, (4T) 



and, K cd = - ./— , 1 3 3 J —1 . (48) 
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which allows the calculation of either phase shift depending on the relative strengths 
of the fields. 

To understand these nonlinearities it is instructive to consider certain limits. One 
of the common limits is when the pump is strong and the single photon fields weak, 
i.e. O2 3> £li^Jn[, fi3 V / n3. In this case the population will be optically pumped 
predominantly into |o) and the N system reverts to two weakly coupled two-state 
transitions. In each case then the phase shift becomes 

i^„ea k = Q 2 J ) . (49) 

which is an ideal cross-Kerr nonlincarity, but is reduced by the Rabi frequency of the 
strong pump. 

Another interesting limit is when Vl\n\ — Q2. Note that this limit will never 
exactly be reached for a coherent state probe due to the uncertainty in m. In this 
case we obtain 

oquai _ _ / ~irnln 3 f A| 

2 7 2 4 V 2 7 ; 

/ A2 



^:r=-A (50) 



(51) 



So there is only a factor of two difference between the schemes. Note that because 
of the probing condition, we should interpret K a b as the desired cross-Kcrr effect to 
realise QND measurement, and K c d represents a self-Kerr effect, although there are 
still cross-Kerr nonlinearities at work as the n\ term has been removed only by the 
special choice of the limit. 

To explore the optimal parameter regime, we define k = fii/f^, which gives 

—fl 3 n 3 J j 



2 7 2""3" a " ^272 J K 2 ni + i' 



^ = -Jo fl 3"8 J k5 , 1 - (53) 



To directly compare these results, in Fig. [5] we show T a b = n 2 ni/ (n 2 ni + l) 2 and 
T c d = K 2 n\/ {n 2 ni + 1) as a function of K 2 n\ for varying n\, which is equivalent to 
the ratio . The T show us the important regimes and there are many features 

of interest here. Firstly we see that T a b is maximised for k 2 — l/^/nT at T a b = 1/4, 
which is the equal case from from Eq. 1501 This demonstrates the result anticipated 
earlier that the largest phase shift is found when the Rabi frequencies of the fields 
in the EIT system are equal. Secondly we note that T ce i is always greater than T a t>, 
and is monotonically increasing. One must realise, however, that the nature of the 
nonlinearity is different for the two probing conditions. By scaling with flini, we are 
implicitly treating the \a) — \b) transition as the probed transition, and \c) — \d) as 
the unknown transition. Hence T a b and T c d in this context refer to scalings of the 
cross-phase and self-phase modulations respectively. Note that T c d asymptotes to 1 
as K 2 ni — > 00. 

The correct way to explore the nonlinear effect and hence QND measurement 
of the unknown field on the probe is to determine the Q-function of the state after 
the interaction. To effect a QND measurement, we require the Q-functions of the 
probe beam with and without a single photon in the channel being monitored to be 
distinguishable. Explicitly, returning to Eq. [46] we take the initial state of the probe 
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Figure 8. Scalings of the cross-Kerr (T a f, - solid lines) and Kerr (T c( j - dotted 
lines) nonlinearities as a function of K 2 n± = for the case of a weak coherent 

pulse applied to the \a) — \ b) transition and the unknown (possibly single photon) 
field applied to the |c) — \d) transition. The maximum of T a (, occurs for the case 
that Q\ = Q2- 



/ I 1 2 \ 00 ni 

exp -£L V ±= exp (iK ab t) \m), (54) 



field to be a coherent state |a), with mean photon number |a| 2 . After interacting with 
this susceptibility, K ab for a period of time i, the probe will be in the state 

\a 

• \ \> 1 

ni=0 

and recall that if a (, is a function of both n\ and 71,3. Note that in general \a'} will not 
be a coherent state following the interaction. We may define the Q-function for the 
state after interacting with the non-linear medium, which is 

Q(/3)H(/3|a')| 2 , (55) 



M 2 +l/3| 2 -2ni , 

...xpl-J-L+lffl 1} /n r X 



^ / |c t | 2 + |/3| 2 -2m 1 \ 

~ expf-^^ M 



exp (iK ab t) x 



V27rmi \ mi 



7) ; 1 



cxp(-z^: h t) (56) 



Note that as required, when ^3 = (i.e. the case of no-photon on the \c) — \d) 
channel), K ab = and the Q-function is simply the same Q-function expected for no 
interaction, i.e. there is no self-Kerr modulation. 

The candidate system for realising a QND measurement that we are considering 
is a monolithic single-mode diamond waveguide (at 637 nm containing NV _ centres, 
which is being pursued using a number of different fabrication strategies e.g. 
[BTH [521 [76] . The energy levels show many possible configurations for achieving the 
four-state system under consideration [161 66] and we will not delve further into these 
schemes, apart from noting that the transition dipolc moments can be achieved and to 
some extent tuned in situ. Maximal coupling requires the smallest width single mode 
waveguides, which for the zero phonon line of NV~ corresponds to a cross section 
of around 200 x 200 nm 2 . To avoid potential cross coupling between centres, we 
require an inter-centre spacing at least around 1 /jm, which for the waveguide under 
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Re[/3] 

Figure 9. Q-functions for various conditions necessary for QND measurements 
using realistic parameters. The three sets of radial distributions correspond to the 
output probe beam when the signal has 0, 1 or 2 photons. Sets of distributions 
at constant radius from the origin are grouped with equal \a). The increase in 
effective resolution with increasing a is demonstrated in this plot, and as discussed 
in the text, there is no apparent self-nonlinear processes. Parameters were chosen 
to correspond closely to expected EIT conditions in low N diamond containing 
NV _ colour centres, and these of reported in Table [T] 



consideration corresponds to an atomic density of TV = 4 x 10 20 m~ 3 . This level 
is dilute but achievable using ion implantation into ultra-low N synthetic diamond 
(e.g. Ref. [77]). Also for NV diamond, we have sr ~ 10 and ix a b ~ 10~ 30 on the 
zero-phonon line transition. 

To determine the bandwidth of the pulse to be measured, we rearrange Eq. [20] 
for an absorption of A — 1%. This implies that the bandwidth of probe and single 
photon field should be / = 8.9 x 10~ 3 r = 740 kHz, or approximately T/100. With this 
bandwidth, we may immediately determine the single photon Rabi frequency, which 
from Eq. 23] is Cl^ = 195MHz. The group velocity depends on the EIT condition and 
also the number of photons in the probe field. 

Under these assumptions, and assuming the reduced inhomogeneous broadening 
for NVD on the zero phonon line that has been observed in low N diamond [TH] [HJ [7H] 
of 7 = 10 GHz, we are now able to fully model the rotation of the probe field in the 
QND measurement. Results from simulation are shown in Fig. [9] which shows the 
Q-function for various values of a and number of photons in the signal beam (field 3). 
A full list of parameters used in the calculations are provided in Table [T] 

Fig.[5]shows families of Q-functions for the state of the probe field after traversing 
the EIT medium. There are three sets of lines. The lowest, horizontal set corresponds 
to the case that the signal field (field 3) had no photons, and represents the 
unperturbed probe field. This highlights the fact that the three-state EIT system does 
not exhibit self-phase modulation. Each point along this set corresponds to increasing 
a from 1 to 40. The upper two curves correspond to the probe field after traversing the 
medium with the signal field in the one photon (middle curve) or two photon (highest 
curve) state respectively, with a increasing to the right as before. QND measurement 
can be inferred whenever the Q-functions of the probe field corresponding to 0, 1 and 2 
signal photons do not overlap. This can be seen clearly in the regime a >~ 15 (fourth 
set of Q-functions). Note that the phase shift is not linear in the regimes we are 
considering, and is greater than would be expected from the simple linear phase shift 
from an ideal cross-Kerr medium. This is likely due to the presence of higher-order 
nonlinearities in the regime Qi — f^, although we have not made a detailed study 
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of them. The distinguishability is maximised at around a ~ 25 for these parameters. 
Note that this is earlier than the maximum expected from T^, which we attribute 
to group velocity dispersion reducing the effective interaction time for certain modes 
of the probe field. A complete study of these processes is beyond the scope of this 
work. A further point to note is that although QND discrimination can be made with 
very high certainty for a >^ 14, there is some noticeable distortion, especially in the 
2-signal photon branch. This indicates that further work is required before strong 
claims of the effectiveness of this scheme for quantum gate operation can be made. 



Table 1. Parameters for NV diamond waveguide system under consideration as 
QND gate 



Waveguide dimensions 


wq x wo x I = 200 nm x 200 nm x 200 /zm 


Atomic concentration 


J\f = 4 x 10 2U m" a 


Dipole moment (ZPL) 


Mah = 10- 30 Cm" 1 


Transition frequency 


A = 637 nm 


Relative permittivity of diamond 


e R = 10 


Homogeneous linewidth 


T = 83 MHz 


Inhomogeneous linewidth 


7 = 10 GHz 


Pulse bandwidth 


/ = 740 kHz 


Semiclassical Rabi frequency 


fi 2 = r/io 


Field ratios 


k = VLi/n 2 = 1/50 


Signal per photon Rabi frequency 


fl 3 = ^VSr/A/Z = 195 MHz 


Detuning scaling 


J = 1/6 


Group velocity 


%.* = "mS" + x ) = 2 - 88 x 104 ($b + x ) ms_i 



The parameters used in the equations used to generate the results shown in Fig. [5] 
are fairly typical of what has already been achieved or will be soon achieved in NVD. 
With variations in the parameters, there will be changes in the resulting phase shifts, 
but the parameters that we chose were not "fine tuned" . 

5. Conclusions 

We have performed investigations of three-state EIT in an inhomogeneously broadened 
sample with the aim of determining the necessary conditions for performing a quantum 
non-demolition measurement in a realistic solid-state medium, diamond containing 
the nitrogen-vacancy colour centre. Our results suggest that even in the presence of 
the relatively large inhomogeneous linewidth of these systems, QND measurements 
are possible using relatively modest extensions of the existing state of the art 
(demonstration of single mode waveguides of order 200 ^m long, and group velocity 
compensated propagation). These conclusions add substantial impetus to the ongoing 
push for diamond based quantum photonics and offer increased support for quantum 
computing based on weak nonlinear interactions. 
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